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ABSTRACT 



Based on simultaneous observations of solar flares in 
hard and soft X-rays we studied several aspects of the 
Neupert effect. About half of 1114 analyzed events 
show a timing behavior consistent with the Neupert 
effect. For these events, a high correlation between 
the soft X-ray peak flux and the hard X-ray fluence 
is obtained, being indicative of electron-beam-driven 
evaporation. However, for about one fourth of the 
events there is strong evidence for an additional heat- 
ing agent other than electron beams. We discuss the 
relevance of these findings with respect to Parker's 
idea of coronal heating by nanoflares. 



1. INTRODUCTION 



The Neupert effect is the name given to the observa- 
tional finding that the rising part of the soft X-ray 
(SXR) light curve often resembles the time integral 
of the hard X-ray (HXR) or microwave emission 
(Neupert, 1968; Dennis & Zarro, 1993). The physical 
relevance of the Neupert effect basically arises from 
the fact that it is interpreted as a causal connection 
between the thermal and nonthermal flare emissions, 
which can be naturally explained within the non- 
thermal thick-target model (Brown, 1971). In this 
model, the flare energy is released primarily in the 
form of nonthermal electrons, and hard X-rays are 
produced via electron-ion bremsstrahlung when the 
electron beams impinge on the lower corona, transi- 
tion region and chromosphere. The model assumes 
that only a small fraction of the electron beam en- 
ergy is lost through radiation; most of the loss is due 
to Coulomb collisions that serve to heat the ambi- 
ent plasma. As a consequence of the rapid energy 
deposition a strong pressure imbalance develops be- 
tween the dense, heated chromosphere and the ten- 
uous corona. The high pressure gradients cause the 
heated plasma to convert into the corona in a pro- 
cess known as chromospheric evaporation (Antonucci 



et al. 1984; Fisher et al. 1985), where it gives rise to 
enhanced SXR emission via thermal bremsstrahlung. 

In this case, the hard X-ray flux is linked to the 
instantaneous rate of energy supplied by electron 
beams, whereas the soft X-ray flux is related to the 
accumulated energy deposited by the same electrons 
up to that time, and we can expect to see the Neu- 
pert effect. Any deviation from the Neupert effect, in 
principle, suggests that the hot SXR emitting plasma 
is not heated exclusively by thermalization of the ac- 
celerated electrons that are responsible for the HXR 
emission. Therefore, investigations of the Neupert 
effect provide insight into the role of nonthermal elec- 
trons for the flare energetics. 



The Neupert effect can be expressed as 
-Pp,sxR = k ■ J-HXR, 



(1) 



with Fp.sxR the SXR peak flux and J^hxr the HXR 
fluence, i.e. the HXR flux integrated over the event 
duration. The coefficient k depends on several fac- 
tors, as, e.g., the magnetic field geometry and the 
viewing angle, and thus may vary from flare to flare 
(Lee et al., 1995). However, if k does not depend 
systematically on the flare intensity, then the SXR 
peak flux and the HXR fluence are linearly related. 



2. DATA SELECTION 



We utilize the SXR data from the Geostationary 
Operational Environmental Satellites (GOES) and 
the HXR data from the Burst and Transient Source 
Experiment (BATSE) aboard the Compton Gamma 
Ray Observatory. The X-ray sensor aboard GOES 
consists of two ion chamber detectors, which provide 
whole-sun X-ray fluxes in the 0.5-4 and 1-8 A wave- 
length bands. BATSE is a whole-sky HXR flux mon- 
itor that, in part, consists of eight large-area detec- 
tors. From each detector there are hard X-ray mea- 
surements in four energy channels, 25-50, 50-100, 
100-300 and >300 keV (Schwartz et al., 1992). 
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Figure 1. Histogram of the difference of the SXR 
maximum and the HXR end time in absolute values 
(top panel) and normalized to the HXR event dura- 
tion (bottom panel). Positive values indicate that the 
maximum of the SXR emission occurs after the end 
of the HXR emission, negative values vice versa. 

For the analysis, the 1-min averaged GOES SXR 
data measured in the 1-8 A channel and the HXR 
data collected in the BATSE Solar Flare Cata- 
log, archived in the Solar Data Analysis Center at 
NASA/Goddard Space Flight Center for the period 
01/1997-06/2000 are used. The peak and total 
count rates are background subtracted for the flux 
below 100 keV. For the SXR events, we used the flux 
just before the flare start for background subtraction. 
To be identified as corresponding events we demand 
that the start time difference between a SXR and 
a HXR event does not exceed 10 min. Overlapping 
events are excluded. Applying these criteria, we ob- 
tained 1114 events that were observed in both hard 
and soft X-rays (for details see Veronig et al., 2002a). 



3. ANALYSIS 

For each event we determined the difference. At, of 
the peak time of the SXR emission, isxR,p, and the 
end time of the HXR emission, ^hxr.e- Further- 
more, the time differences were normalized to the 
duration D of the respective HXR event, i.e. 

A , _ ^ _ ^SXR.P - tHXR.E 

^tnorm — „ — „ . (.^J 
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Figure 2. Scatter plot of the SXR peak flux versus 
the HXR fluence. The regions that lie beyond the es- 
timated HXR fluence and SXR peak flux thresholds 
are grey shaded. Moreover, we have indicated the 
sign of the time difference between SXR peak and 
HXR end for each single event: symbols de- 

note events with positive, "— " symbols events with 
negative, "o" symbols events with zero time differ- 
ence (i.e. the SXR peak and HXR end take place 
within Imin). The straight line indicates a line of 
constant k, i.e. i^p,sxR = 2 • 10^^^ • J-hxr- 



Figure ^ shows the histogram of the absolute and 
normalized time differences. Both representations 
of the SXR-HXR time difference have its mode 
at zero. 49% of the events lie within the range 
\M\ < 1 min, and 65% within \M\ < 2 min. For 
the normalized time differences, we obtain that 44% 
lie within |Atnorm| < 0.5 HXR units, and 59% within 
|At,ioim| < 1 HXR unit. This outcome suggests that 
certainly a considerable part of the events coincides 
well with the expectations from the Neupert effect 
regarding the relative timing. 

Figure ^ shows the scatter plot of the SXR peak 
flux versus the HXR fluence for the complete sample, 
clearly revealing an increase of fp^sxR with increas- 
ing J?^HXR- It can also be inferred from the figure that 
the slope is not constant over the whole range but 
that it is larger for large HXR fluences than for small 
ones. For very large fluences, the slope approaches 
the value of 1, indicative of a linear relation between 
the SXR peak flux and HXR fluence. We stress that 
the slope at small fluences might be affected by miss- 
ing events with small SXR peak fluxes (due to selec- 
tion effects), and thus appear flatter than it is in fact. 
Figure || reveals an interdependence between the im- 
portance of an event and the sign of the time differ- 
ence. Basically all large flares belong to the group 
of events with At < 0, i.e. the SXR peak occurs be- 
fore the HXR end. Moreover, the flares with Ai < 
reveal a strong tendency to be of long duration. 

We obtain a high cross-correlation coefficient for 
the SXR peak flux and HXR fluence relationship, 
r = 0.71. This coefficient is higher than those for 
the SXR peak flux and HXR peak flux, r = 0.57. 
This indicates that the correlation is primarily due 
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to the HXR fluence - SXR peak flux relationship, 
as predicted from the Neupert effect, and not, e.g., 
due to the fact that flares with large HXR peak fluxes 
also tend to have intense SXR counterparts. Further- 
more, it is important to note that the HXR fluence - 
SXR peak flux correlation is higher for the events 
with negative time differences, r = 0.82, than for the 
events with positive time differences, r — 0.54. 

On the basis of the relative timing of the SXR 
peak and the HXR end, we extracted two subsets 
of events. The events of set 1 are roughly consistent 
with the timing expectations of the Neupert effect, 
and the events of set 2 are inconsistent with it. The 
two sets are defined by the following conditions: 

Set 1 : (|At| < 1 min) OR (|Ainorm| < 0.5 unit) , 
Set 2 : (|Ai| > 2 min) AND (|At„orin| > 1-0 unit) . 

The applied conditions represent a combination of 
absolute and normalized time differences in order 
to avoid as much as possible any a priori interde- 
pendence with the flare duration and/or flare inten- 
sity. Out of the 1114 corresponding HXR and SXR 
flares, 485 (44%) events fulfilled the timing criterion 
of set 1; 270 events (24%) belong to set 2; 359 events 
(32%) are neither attributed to set 1 nor to set 2. 

In Figure H, we plot the SXR peak flux versus the 
HXR fluence separately for set 1 and set 2. The 
figure reveals that the two sets have very different 
characteristics besides the different timing behavior. 
Set 1 contains many more large events and shows 
a steeper increase of i^p,sxR with increasing J^hxr 
than set 2. Moreover, set 1 contains more events 
with negative than positive time difference, whereas 
almost all events of set 2 are characterized by At > 0, 
i.e. increasing SXR emission beyond the end of the 
hard X-rays. On average, for small fluences the 
events belonging to set 2 have a larger SXR peak 
flux at a given HXR fluence than those of set 1, in- 
dicating an "excess" of SXR emission with respect 
to set 1. The cross-correlation coefficients derived 
separately for the subsets reveal that the correlation 
among the SXR peak flux and HXR fluence is much 
more pronounced for the events of set 1, r = 0.78, 
than those of set 2, r = 0.41. 



4. DISCUSSION AND CONCLUSIONS 

24% of the events have Ai < 0, i.e. the SXR max- 
imum occurs before the end of the HXR emission. 
These events are preferentially of long duration. Li 
et al. (1993) have calculated time proflles of soft and 
hard X-ray emission from a thick-target electron- 
heated model, finding that, in general, the time 
derivative of the SXR emission corresponds to the 
time profile of the HXR emission, as stated by the 
Neupert effect. However, for gradual events they ob- 
tained that this relationship breaks down during the 
decay phase of the HXR event, in that the maximum 
of the SXR emission occurs before the end of the 
HXR emission. This phenomenon can be explained 
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Figure 3. Scatter plot of the SXR peak flux versus the 
HXR fluence separately plotted for set 1 (top panel) 
and set 2 (bottom panel). The straight lines are lines 
of constant fc (= 2 • lO'^^)^ 

by the fact that the SXR emission starts to decrease 
if the evaporation-driven energy supply cannot over- 
come the instantaneous cooling of the hot plasma, 
which is likely to happen in gradual flares. Con- 
sidering our observational flndings together with the 
results from simulations by Li et al. (1993), presum- 
ably most of the events with Ai < are consistent 
with the electron-beam-driven evaporation model. 
In particular, the high correlation between i^p,sxR 
and J^HXR, ^ ~ 0.8, supports such interpretation. 

56% of the events have At > 0; these events are 
preferentially of short and weak HXR emission. In 
principle, the fact that the SXR emission is still in- 
creasing although the HXR emission, i.e. the elec- 
tron input, has already stopped indicates that an 
additional agent besides the HXR emitting electrons 
is contributing to the energy input and prolonging 
the heating and/or evaporation. However, McTier- 
nan et al. (1999) have shown that the SXR time 
profile depends on the temperature response of the 
used detector: An increase of the SXR emission 
of low-temperature flare plasma after the HXR end 
may arise due to cooling of high-temperature plasma. 
Thus, we cannot attribute all events with At > as 
inconsistent with the electron-beam-driven chromo- 
spheric evaporation model. Instead, we consider as 
inconsistent only flares which show strong deviations 
from At = 0, i.e. the events belonging to set 2. 
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Figure 4- Flare frequency distributions as function 
of the SXR peak flux (top panel) and HXR fluence 
(bottom panel). The straight line indicates the linear 
fit in log-log space, characterized by a slope —a. 



This means that for at least one fourth of the ana- 
lyzed events an additional heating agent besides non- 
thermal electrons is suggested. A probable scenario 
is that energy is transported from the primary en- 
ergy release site via thermal conduction fronts, which 
initiate chromospheric evaporation but do not pro- 
duce hard X-rays. The finding that for a consider- 
able fraction of flares, preferentially weak ones, an 
additional heating agent other than electron beams 
is suggested, is not only relevant for the flare ener- 
getics but also for Parker's idea of coronal heating 
by nanoflares (Parker, 1988). 

Hudson (1991) pointed out that if the corona is 
heated by flare-like events of different sizes, then the 
flare energy distribution must have a power-law slope 
a > 2. If the SXR flux does not vary systematically 
with temperature and density, then the SXR peak 
flux is linearly related to the maximum thermal en- 
ergy of the flare plasma (see Lee et al., 1995; Veronig 
et al., 2002a). On the other hand, HXR fluence dis- 
tributions can be considered as representative for the 
energy contained in nonthermal electrons. In Fig- 
ure El we show the SXR peak flux and the HXR 
fluence distributions derived from 1114 correspond- 
ing SXR/HXR flares, flnding a = 1.95 ±0.13 and 
a = 1.58 ±0.13, respectively. The discrepancy be- 
tween the slopes of the HXR fluence and the SXR 
peak flux distributions was already pointed out and 
discussed in Lee et al. (1993, 1995) and Veronig et 
al. (2002b). The present analysis provides an expla- 



nation for this difference in power-law slopes: The 
relationship between the SXR peak flux and the 
HXR fluence is not linear, whereby the deviations 
from a linear correlation are strongest for weak flares 
(cf. Figure ||) . 

The soft X-ray flare emission increases due to energy 
supply by electron beams as well as any other heating 
agent, whereas the hard X-ray emission contains only 
information on the energy provided by electrons. To- 
gether with our finding that particularly in weak 
events an additional heating agent besides electron 
beams is suggested, this strongly suggests that soft 
X-ray peak flux distributions are a more meaning- 
ful indicator of flare energy distributions than hard 
X-ray fluence distributions. Furthermore, we have 
shown that weak flares have different characteristics 
than large flares, in the sense that electrons are less 
important for their energetics. Therefore, it is possi- 
ble that the power-law slope of nanoflare frequency 
distributions differs from that derived for observed 
flares, which is close to the critical value of 2. In this 
respect it should be worthwhile to investigate flare 
frequency distributions of SXR flares without HXR 
counterparts, i.e. without detectable particle acceler- 
ation, since these flares possibly provide the link to 
the smallest flare-like energy release events. 
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